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INTRODUCTION 


It is common to express the bending fatigue life of a material in 
terms of the nominal elastic stress or strain calculated by assuming 
stress is proportional to strain. Yet, in the low-cycle fatigue range, 
and in fact over a wide range of life of engineering interest, plastic 
flow is involved, and stress is not proportional to strain. Thus con- 
siderable discrepancy exists between bending fatigue and axial fatigue 
when the latter involves expression in terms of true stresses and strains. 
While several factors contribute to this discrepancy, as discussed in 
Ref. 1, among them the presence of strain gradient, a volumetric effect, 
cyclic strain hardening and softening effects, and crack propagation 
differences after the surface stress has been initiated, the major effect 
is due to the fact that the nominal stress and strains on a specimen in 
bending differ considerably from the true values due to plastic 
deformation. 

In Ref. 1, it was demonstrated that flexural fatigue is consider- 
ably different from axial fatigue, but that the two could be brought in‘o 
coincidence when based on true surface stresses. Fig. 1 shows some of 
the results of the reference. The solid curve and circular data points 
show the axial fatigue characteristic of 4130 steel. In this case the 
stresses are true values because in axial loading it is simple to deter- 
mine the stress from a knowledge of load and cross-sectional area. The 
dashed curve and square data points refer to the beam bending 
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characteristics. Here the stress is calculated from the conventional 
Me 

S = T assumptions, so that the stresses are the fictitious values 
determined on the basis that no plastic flow occurs. It is seen that 
for this material the nominal bending stress can be 50% higher at the 
100 cycle life level, and that significant difference persists between 
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the two curves until well above 10 cycles to failure. Obviously, then, 
it is important to reconcile the differences over a wide range of life 
levels of practical engineering interest. It wfs also demonstrated in 
Ref. 1, that if the true stress is calculated in the case of the bending 
specimen the li*e comparison with the axially loaded comparison is much 
more favorable. But, because of the difficulty of performing the inte- 
grations involved in the bending case no closed-form solution was 
presented. Each material and geometry had to be treated numerically. 

It is the purpose of this report to extend the work of Ref. 1 and 
to obtain closed-form generalized relationships. The approach is first 
tc study rectangular rather than circular sections. Not only are 
rectangular sections of great importance in engineering applications, 
but they result in relations that can be integrated in closed-form. 

The resulting form of the solutions then provides a model for use in 
the case of circular sections. By assuming an analogous model, 
involving several undetermined constants, and then determining the 
constants so as to fit a range of known solutions obtained by numerical 
integration, it has been possible to obtain very accurate, though not 
exact, closed-form solutions for the case of the circular section. 

The procedure and results will be described in this report. 
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In addition to the development of plasticity as a reason for 
difference between axial and bending fatigue, another factor relating 
to volume under stress also introduces differences in the lives of the 
two loading modes. In pure axial loading the entire cross section of 
the specimen is subjected to the loading stress, whereas in pure bend- 
ing only the surface fibers are subjected to the maximum stress, wriile 
other regions of the cross section are subjected to lower stresses. 

The effect of this factor is expected to be important in the high cycle 
fatigue range, rather than the low cycle fatigue range which is of major 
interest in this report. However, this subject is also briefly 
di scussed. 


BASIC EQUATIONS USED FOR LIFE RELATIONSHIPS 

The Conventional Life Relationship 

The basic model for axial fatigue used in this report is shown in 
Fig. 2. Although the model was first proposed by Manson (Ref. 2) 
using a different notation, the notation used in Fig. 2 is the one that 

A, 

has been commonly adopted ( Ref. 3). The strain amplitude, y is 
plotted on the vertical axis, and the number of reversals 2N^ on the 
horizontal axis, using logarithmic scales on both axes. If a material 
is axially cycled at constant strain range about a zero mean value, or 
if the plastic strainrange is maintained constant, it is generally 
found that in the early cycles the stress (or, equivalently, the elastic 
strain) may change considerably, but eventually a saturation value is 
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achieved. Usually this saturation value develops early and persists 
for the major part of the life of the specimen. Plots such as shown in 
Fig. 2, in which both the plastic strain amplitude, and the elastic 
strain amplitude, are plotted against reversals to failure result in 
straight lines PQ and LM, the equations of which are, respectively, 

Y = e f( 2 N f) c (1) 

Ae . al . 

—r’Ht’T < 2 N f> (2) 

where Ae is the plastic strainrange, Ac is the elastic strainrange, 

P , oi 

and Ac the saturation stress range, e^. and 3re the intercepts 

at 2N^ = 1 of the plastic and elastic lines, E is the elastic modulus, and 

c and b are the slopes of the two lines as shown in Fig. 2. 

The equation for the total strain amplitude, which is the sum of the 

plastic and elastic components, is thus 

Q l 

Y = e}(2N f ) c + (2N f ) b (3) 

In Ref. 4 we introduced the term "transition life" to designate the 
life at the intersection point T of the elastic and plastic lines. In 

Fig. 2 this life is denoted Ny and the number of reversals is 2Ny. 
Correspondingly, the strainrange at this point is Ae^ and the strain 
amplitude is Ae^* By equating the strains in Eqs. (1) and (2) the 
values of Ny and A tj can easily be determined. 
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N 


T 



1 

b-c 



•v 


(4) 


(5) 


where n = ^ . It may be noted that n is the reciprocal of the strain 

rr, * 

hardening exponent m' in the relation Ao = Constant (Ac p ) . From 
Eqs. (4) and (5) it can readily be shown that the basic life model 
Eq. (3) can be rewritten [5] 



The Inverted Life Relation 

For convenience in some of the numerical integration analyses we 
have also used an alternative form of the life relation discussed in 
Ref. 6. Basically, life is expressed directly in terms of strain, 
instead of using Eq. (6) which expresses strain in terms of life 




+ R 


(7) 


where 
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z = exp 


P ln 2 R £ +Q In R £ + ln|-.889c (|) 36 j 


( 8 ) 

(9) 

( 10 ) 


P = -.001277 (|) + .03893 (|) - .0927 

c 2 c 

Q = +.004176 (|) - .135 (^) + .2309 

R £ = — = ratio of strain required to produce life Ny to 

the transition strain (11) 


For any material of known b, c, Ny and Aty, all the constants 
involved in the eq* ’f ions are readily calculated, the life Ny is ex- 
pressed directly in terms of those constants and strainrange. 

Similarly, stress range An can be written directly in terms of 
applied strainrange Ac 


A o = EACy 


R z/b 



1 b/z 


( 12 ) 


The Universal Slopes Life Equation : In cases where material 

constants b, c, e'y and o'y are not accurately known from experiment, 
it may be desired to estimate the life relation on the basis of common- 
ly available material properties. In Ref. 1 it was demonstrated that 
as a first approximation c could be taken as a universal value of -0.6, 
b as a universal value of -1.12, and that e'y could be estimated from 
the ductility, while o'y was best estimated from ultimate tensile 
strength. The resulting "universal-slopes" equation, as expressed in 
Ref. (1). is*, 




3 * 5o u t. -0.12 
E N f 


(13) 


★ 

While this equation is generally credited only to Manson because 
it first appeared in Ref. 1, the equal contribution of M. H. Hirschberg 
in its development is acknowledged. 
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where 

D = ductility = -In (1-RA) wherein RA is reduction of area in v % 

tensile test 

o = ultimate tensile strength, and E is elastic modulus in 
u consistent units. 


Expressed with the SAE notation of Eq. (3), the equation becomes 


4r = 0. 7578D°' 6 (2N f )" 0 ' 6 + 1.902 ^ (2N f )" 0,12 

(14) 

That is, c' f = 0.7578D 0,6 and a‘ f = 1.902° u 

(15) 

Also, using Eqs. (4) and (5), 


Ae-j. = 4.789 (^H) 1 ' 25 D"' 15 

(16) 

*^.083 

N t = .0735 (-g-) D 1 '*’ 

(17) 


so that Eq. (6) becomes 



where Acy and N-. are given by Eqs. (16) and (17), respectively. 

Similarly, substituting c = -0.6, b = -0.12 into Eqs. (9) and 
(10), we obtain P = .07, Q = -.3397 

z = exp[.07 ln 2 R e -.3397 lnR £ -1.208] 



N f = N t [R l z/c + R e z/b ] 1/z (19) 

where Atj, Ny are given by Eqs. (16) and (17), respectively 
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FLEXURAL BENDING OF RECTANGULAR CROSS-SECTIONS 
We first treat the bending of rectangular cross-sections because 
the resulting equations can be integrated in closed form. As discussed 
in Ref. 1, it is more convenient to solve this type of problem in an 
inverse fashion compared to conventional treatment. Rather than start- 
ing with a known bending moment, and determining the resulting surface 
strain--which would involve solution of non-linear equations- -it has been 
found better to select a surface strain and calculate the oending moment 
required to produce this strain. From the selected surface strain we 
establish the fatigue life, and from the required bending moment the 
nominal elastic stress. Thus we can establish the relation between nominal 
elastic stress and fatigue life. 

Fig. 3 shows the notation used in the analysis. We assume a sur- 
face strain corresponding to a life N<. when the full bending moment M 
is applied. Since in bending it is common to assume that plane sections 
remain plane, strain varies linearly across the thickness, so that the 
strain at a distance y from the neutral axis (center of bar) is 


e 


y 



( 20 ) 


Substituting for t in terms of life N $ associated with the surface 
strain, from Eq. (6)* 


c 


y 



( 21 ) 


♦Because’ "of the symmetry or loading for this ose of completely 
reversed bending, we omit the A's, and consider only surface strains 
at the maximum loading condition. 
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Now if Ny is the life associated with the strain at y, then by Eq. (6) 


N \ b IH \ c 
-X- 1 + JM 

N T \ T / 


and the stress associated with the strain is, as derived from Eqs. 
( 2 ), (4) and (5), 


ECt K 


The increment of bending moment contributed by a strip w dy at toe 
distance y from the neutral axis becomes 


= y • Oy w dy - E e T ( j * w y dy 


and the total bending moment, taking into account the symmetrical contri- 
bution to bending moment of the section below the neutral axis is 

f y=>1 b 

M = JdM = 2Ewe T ||^ j y dy (25) 


Here N y is a function of y according to Eqs. (21) and (22). In order 

to carry out the integration it is best to establish dy in terms of 
N N 

fit , rather than expressing j?- in terms of y, which would result in a 
™T n T 

more difficult integration. Thus from Eq. (21) and (22), 

(%)“•(%)' 

y • •> -h- (b re <»: 

i Mi 
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<30 



b 

3S 



b+c 

2b+c 



2b+c 



( 30 ) 


In Eqs. (28), (29) and (30) the lower limit for ^ is taken at 
h = 0, where the strain is zero, and therefore the life is infinite. 



r 


Since all 
the value 


N 

the exponents of if- in the integrated expression are negative. 

i 

of all terms is zero at the lower limit, and M is evaluated a 1 



Mh 


The nominal surface stress S. .. is f 

bending I 


12Mh 


1 2M 


since the 


w(2h) 8wh‘ 

total height of the rectangle is 2h and the moment of inertia of a 

1 3 

rectangular section is base x (height) . Thus, substituting for K 

it - \2b 


from Eq. (31), and dividing both numerator and denominator by 
the resulting equation for S ben(J ^ n g is obtained as follows: 







4 ”- - 3 s correction factor * 6 ( 34 ) 

^axial 


(c-b) 2(c-h) 



6 * 


(35) 



where 


f _ 3(b+c) 1 + n 

r l " _ 2Tj+c~ " W3) + (l/3)n where n = c/b (36) 

f _ _3c_ _ 3n 

2 b+2c l+2n (37) 

FLEXURAL BENDING OF CIRCULAR CROSS-SECTIONS 
We calculate the bending moment required to produce a selected 
surface strain. From the selected surface strain fatigue life can be 
established and similarly nominal stresses from the required bending 
moment. Thus we establish the relation between nominal elastic stresses 
and fatigue life. 

Figure 4 shows the notation used in the analysis. We assume a 
surface strain c corresponding to a life N^ when the full bending 
moment H is applied. Strain varies linearly across the section, so that 
the strain c at a distance y from the neutral axis is 

y 

t = * e 

y R s 

y = R sine , dy = R cosede (38) 

e y = e $ sine (39) 


and the stress associated with the strain 

T / K 


can be found using Eq. 
//I b/z 


( 12 ), 
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where z is given by Eqs. (8) to (11). 

The increment of bending moment contributed by a strip at the 
distance y from the neutral axis becomes 


dM = y • o • (width of the strip) *dy 
Substituting for y and dy from Eq. (38) and from Eq. (40), 


dM = R sin0 • Ee-, 


e s sine ^ sine^^ 


z/b 


2Rcose Rcosedi. 
(41) 


The total bending moment, taking into account the symmetrical contribu- 
tion to bending moment of the section below the neutral axis is 



where z is a function of (c/b) and U s /e-j.) given by Eqs. (8)-(ll). 

It is clear from Eq. (42) that the integral does not lend itself to 
simple closed form solution. However, for any given material and a given 
surface strain, the integral can be computed numerically and the moment 
can be calculated. The nominal surface stress is expressed as: 


X 
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c _ M_R M R 64 _ 64_M 

bendi ng " ' I ‘ ^4 ' ^3 


The ratio of nominal surface stress S to the stress that would be 
required to produce a life N<. in axial fatigue, can be written analogous 
to Eq. (35) by assuming the denominator and the exponent (c-b) remain 
the same for both rectangular and circular cross-sections. 



where f-| and f,, are for the circular cross-section. The ratio 6 is calcu- 
lated numerically using Eqs. (43) and (42) tor different values and 
combinations of parameters involved as follows: 


N 

N 


S 

T 


10 ' 3 , 10 ‘ 2 ,.. . , 10 2 , 10 3 


c = -0.3, 0.4,. .. ,-0. 9,-1. 0 

b = -. 05,- . 06, ...,-. 1 9,-. 20 


Equation (44) can be written as 

^N, 


1 + 21 r 

n t/ 


\(c-b) / 1 . \2(c-b) 

' - 4 ) 

(c-b) /N \2(c-b) 


No\^ C ” b ^ 

= 1+f llN f I + f 2N^ 


(45a) 




% 


Ic-b) 


- 1 


x (c-b) 


= f l + f 2 l N, 


(45b) 
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( 45 ) 


If the function G is plotted against the function H for a given 

N S 

(c/b), i.e. n, but for all combinations of c,b, and , a straight 

T 

line results. Thus the slope f£ and intercept f-j are functions of only n. 

The plot of G vs. H for different n values are shown in Fig. 5. For clan 1 -' 
the numerical values of functions G and H for different combinations of 
c, b and N^/Nj are shown in Fig. 5 only for the case of n = 5. The above 
procedure is repeated for several materials of various n values and the 
parameters fj and f£ are calculated for each material with a specific n 
They are shown in Fig. 6. 

Since the parameters f^ and f^ are found to be functions of only n, 
an attempt was made to determine these functions whose form is analogous 
to ones shown for rectangular cross-section in Eqs. (36) and (37). 

The functions were determined by least squares analysis to be 


1 + .8432 n 
1 " .6849+. 255 n 


(46) 


2.6188 n 

2 " 1+1. 5411 n 


(47) 


The curves represented by the above functions are shown by solid lines 
in Figure 6. 

For any given material the functions f^ and f^ can be calculated using 
the reciprocal of the strain hardening exponent, n, and Eqs. (46) and (47). 
The correction factor 6 and axial stress S axial can be calculated for any 
given life using Eqs. (44) and (33), respectively. Thus we can plot nomi- 
nal bending stress Sb enC |i n g against life. The plots are shown in Fig. 7 
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for several materials of engineering interest. The properties of these 
materials, obtained from Ref. 7, are listed in Table II. 

If the fatigue properties of the material are not known or cannot 
be obtained easily. Universal Slopes solutions can be used with the 
properties obtained from a tensile test. Universal plastic and elastic 
exponents of -0.6 and -0.12 are used, respectively. Transition strain 
V:. T and transition life N y can be determined from Eqs. (16) and (17), 
respectively. 


DISCUSSION 

The approach presented in this report takes into account only the 
effect of plastic flow in explaining the difference between axial and 
bending fatigue lives; hence, it is valid in the low cyclic life range. 
According to the analysis, the S axial and S bending curves become coinci- 
dent for most engineering materials at lives greater than 10^ cycles 
since the strains are purely elastic at thoselife levels and the Mc/I 
assumption is valid. However, bending fatigue life could still be 
higher than at axial fatigue for the same nominal stress due to volumetric 
effect which is not taken into account in the current approach. This is 
observed by Esin in steels and he has used micro-plastic strain energy 
criterion to take this factor into account (Ref. 8). This is also 
diseased in some detail in connection with risk of-rupture criterion 
In Ref. 9. Basically, in bending fatigue lower volume of material is 
under high strain at outer fibers as compared to the whole volume in 
the axial case. Thus, the smaller the volume under maximum stress. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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the lower will be the probability of the presence of a flaw at a point 
of high stress, hence the higher will be the life. 

At high cyclic life range, the rotating bending will give lower life 
as compared to alternating bending for the same nominal stress which can 
e explained as follows. In alternating bending the regions A and B in 
Fig. 4 will be subjected to maximum surface strain, whereas in rotating 
bending all the material near the surface of the cylinder will be subject 
to maximum strain at different times; for example, the region C will take 
the position of A or B while rotating, hence will be subjected to maximum 
strain. Thus, more volume is subjected to high stress in rotating bending 
hence greater will be the probability of the presence of a flaw at a point 
of high stress and lower will be the life. 

It should be noted that results obtained in this report are valid 
for alternating flexural bending. If the material is subjected to rotating 
bending, the cross-section deflects about a different axis, which makes an 
angle with the loading axis ( Ref. 1). Once the angle between two axes i: 
determined, the complete stress and strain distribution can be computed 
from which the bending moment for an assumed value of maximum surface strain 
can be established. In Ref. 1 i* was observed that although the analysis 
for rotating bending is considerably more complicated than static flexural 
bending, the resulting strain distributions for the two cases differed 
little when considered in their respective bending planes. It was shown 
that the experimental results for rotating bending agreed well with the 
predictions made for static flexural bending. Hence it will be assumed here 
that the equations for flexural bending, Eq. (44) together with the f^ and f 2 
values given in Eqs. (46) and (47), are also valid for rotating bending in 
the low cycle fatigue range. Further study may be desirable to verify this 
assumption. 



CONCLUSION 


Closed-form expressions were obtained for nominal bending stress 
in terms of axial fatigue stress based on true surface stresses. They 
are shown for rectangular and circular cross-sections in Table I. Using 
only axial fatigue data on a given material and formulas given in Table I, 
nominal bending stress-life curve can be established. Though the study 
is hunted to only rectangular and circular cross-sections, the approach 
and the procedure described in this report can be applied to more compli- 
cated geometries to obtain the functional forms for parameters f^ and f ^ 
which are functions of only strain hardening exponent of the material. 

The approach takes into account only the effect due to plastic flow which 
has dominating effect in the low cyclic life range. However, in the 
high cyclic life range other factors could be more dominating like the 
volumetric effect which should be taken into account in reconciling the 
difference between axial and bending fatigue lives. Frture study is 
necessary to incorporate the volumetric effect in the analyses presented 
in this report. 




G»' pOGK ^ 
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TABLE I 

n - c/b ratio of plastic to elastic exponents 

1 

k-j = transition strain range = 2c^'^ 1 

where = fatigue ductility coefficient 

o f ' r fatigue strength coefficient 

E = modulus of elasticity 



) n 


S 


bending 






axial 


N 


where f, and f„ are functions of n and depend on the geometry of the 


cross-section. 

For Rectangular Section 

* 1 + n 

f l = (2/3) + (l/3)n 

r . 3n 
f 2 ' TTTH 


For Circular Cross-section 

* 1 + .8432 n 

T 1 " .6849+ .255 n 

, 2.6188 n 

f 2 = T+ T75"4TTn 


For any given surface strain c $ , surface life N c can be found from Fig. 


'‘axial and S 


bending 


can be computed using the above formulas. 


Cvl 
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Table II (Continued) 
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Fig. 1. Fatigue data under axial loading and in rotating bending 
for 4130 steel from Ref. 1. 



strain amplitude 





fig. 4. Circular cross-sectional bar in flexural bending. 




O SAE 1015 - 80 BHN 
-O- VAN 80 - 225 BHN 
<> RQC 100 - 298 BHN 
D SAE 1045 - 500 BHN 
-D- SAE 4142 - 380 BHN 
£ AISI 4340 - 409 BHN 
□ AISI 52100 - 518 BHN 
O- H -11 - 660 BHN 

0 SAE 9262 - 280 BHN 
A AISI 304 - 327 BHN 

-A- 18% Ni Maraging 460 BHN 

1 AM 350 - 496 BHN _/ 


^ = Eq. (46) 


= Eq. (47) 


0 2024 T 351 A1 
-0- Manten - 150 BHN 
$ AISI 4130 - 258 BHN 
■ SAE 950X - 150 BHN 

A SAE 4142 - 450 BHN 


Fig. 6. The parameters fi and f- are shown for materials with 
different n values. The solid lines are curve-fits 
given by Eqs. (46) and (47), respectively. 
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Fig, 7 (Continued) 





‘binding 



Fig. 7 (Continued) 
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Fig. 7 (Continual) 




